Ab-initio band structure of doped graphene 
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We present an ab-initio study of the graphene quasi-particle band structure as function of the 
doping in GW approximation. We show that the LDA Fermi velocity is substantially renormalizcd 
and this renormalization rapidly decreases as function of the doping. We found, in agreement with 
previous papers, that close to the Dirac point the linear dispersion of the bands is broken but this 
behavior disappears with a small doping. We discuss our results in the light of recent experiments 
on graphene and intercalate graphite. 
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I. INTRODUCTION 



II. METHODS 



Since its discovery [l|, 0, H[ graphene raised a great 
interest in the scientific community for the possibility 
to use it as new material with tunable properties |4( in 
electronic devices and to probe open quantum mechan- 
ics questions. Electron-electron interaction in graphene, 
and related compounds, are expected to play an impor- 
tant role due to its low dimensionality. Lately the quasi- 
particle dynamics of graphene and graphite has been ad- 
dressed by high-resolution angle resolved photoemission 
spectroscopy ARPES 0, 0, Q experiments. It has been 
found that conical bands are distorted due to many-body 
interactions, which renormalizes the Fermi velocity vf- 
The self-energy corrections to the bare band structure are 
crucial for determine transport, optical properties and 
other properties of graphene based materials. 
Moreover recent experiments on doped graphene Q re- 
port variation of the Fermi velocity and Fermi surface 
topology induced by doping. Minor changes are definitely 
due to the small variation of the lattice constant Q with 
the doping but mainly it is expected that larger differ- 
ences are due to correlation effects in doped graphene 
fiol . ITl| . In this paper we study the quasi-particle band 
structure of extrinsic graphene (doped graphene) as func- 
tion of the doping, gated or electrostatic. Correlation ef- 
fects are evaluated using ab-initio many-body self-energy 
calculated in GW approximation, where the polariza- 
tion needed to evaluate W is calculated within the Ran- 
dom Phase Approximation. The validity of the RPA in 
the calculation of graphene self-energy has been recently 
questioned [e| but it has successively shown to be a valid 
and controlled approximation in doped graphene for any 
sufficient small doping [Io| . Therefore in this paper we 
always introduce a small broadening in the Fermi distri- 
bution of the electrons even in the undoped case in order 
to guarantee the validity of our approximations and be 
closer to the experimental results where a finite doping 
is always present, due to impurities or substrate interac- 
tion. 



In order to simulate isolated graphene we used a slab- 
geometry, i.e., bulk geometry with large distance between 
the layers, and an elementary 2D graphene cell with 2 
atoms with a lattice constant a = 2.46 Angstrom. DFT- 
LDA calculations were performed with the plane wave- 
code ABINIT [l3|, using a cutoff of 60Ry for the wave- 
function, a 42x42a; 1 Monkhosrt-Pack grid for the density 
and Troullier Martins pseudo-potentials [l4[ for the car- 
bon atom. The quasi-particle band structure were suc- 
cessively obtained from the self energy T,(E) evaluated 
in the GW approximation 0, [l6|. GW calculations 
were performed with the code Yambo (l7j starting from 
DFT-LDA wave- functions. Because we are not interested 
in the total energy but only in spectral properties, we 
didn't perform a self-consisted calculation but stop this 
approximation to the first order £ = Go Wo [HI. The 
quasi-particle band structure is then obtained as: 



(1) 



where AE = £ — V xc , V xc is the LDA exchange- 
correlation functional, e n k are the Kohn-Sham eigen- 
values and Z is the renormalization factor Z n k = [1 — 
dReA'S/dw]" 1 . The screened electron-electron interac- 
tion W have been calculated within random phase ap- 
proximation (RPA) in term of the dielectric constants 
€G_Gi(q,uj) using a plasmon-pole model Bands up 

to 60eV in energy have been employed both for W and 
for the expansion of the Green Function G, that corre- 
spond roughly to 50 bands in the case of 20 a.u. distance 
between the graphene sheets. Convergence in the size 
of the matrix eG,G'(q,u) has been carefully checked. In 
all the GW calculations a small Fermi-Dirac broadening 
0, 002i?y was used. The dependence of the GW results 
from the k-point sampling and the distance between the 
sheets is discussed in the Appendix A. 
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LDA 


GW 


KM 


0.802(1) 


1.00(5) 


KT 


0.860(2) 


1.15(5) 



TABLE I: Fermi velocity of graphene extrapolated for infinite 
distance between the slabs along the direction KT and KM . 



III. QUASI PARTICLE BAND STRUCTURE OF 
GRAPHENE 



We report the quasi-particle band structure of 
graphene in fig. [1] for the tt and it* bands along the 
symmetry lines TKMT in LDA and GW approxima- 
tion. It's clear that the bandwidth is increased in GW 
and consequently the Fermi velocity vf is renormalized. 
We report vp along the direction KT and KM in the 
table I, extrapolated for an isolated slab of graphene (see 
Appendix A). 
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FIG. 2: Band Structure of Graphene in LDA and GW ap- 
proximation close to the Jfpoint . (a) Undoped graphene, 
(b) doped graphene with 0.05 electrons in excess in the prim- 
itive cell. 
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FIG. 1: Band Structure of Graphene in LDA and GW ap- 
proximation, using 20 a.u. slab geometry and 36a;36a;l kpoint 
sampling. 



Very close to K we found a non-linear behavior of 
the 7r bands similar to the one obtained by Trevisanutto 
et al. ref. [2(|, see fig. [2] panel a). However we notice 
that the slope of the bands we obtained close to K is 
slightly different and this can be due to the plasmon-pole 
approximation employed in the calculations, see also ref. 
[2Cj | . We then studied the anisotropy of the bands close 
to K in GW and LDA. We evaluate the Fermi velocity in 
different directions and found that the GW corrections 
increase the anisotropy respect to the LDA in undoped 
graphene. The results for vf in different directions 
extrapolated for isolated graphene are reported in table 
III 



IV. THE EFFECT OF DOPING 

Doped graphene is simulated changing the number of 
electrons in the unit cell and then compensating the neg- 
ative charge with an uniform positive background. When 
we start to dope graphene the system evolves from a 
semi-mental to a real metal, and the effective interac- 
tion felt by the electrons starts to be weaker because of 
the stronger screening of the Coulomb potential. This 
causes a smaller renormalization of the band structure 
and a consequent decrease of the Fermi velocity. In fig. 
|3] we report the gap at the M point of graphene versus 
the electrostatic doping, that is represented by An the 
excess of electron in the unit cell, divided for the LDA 
gap at zero doping. As one can see while in LDA the gap 
slightly varies with the doping, on the other hand the 
GW renormalization is very strong for intrinsic graphene 
and decreases rapidly as function of the doping. 
Then we analyzed the band structure close to the Dirac 
point and we found that the kink rapidly disappears with 
doping. We report in the figUl the two extreme cases the 
intrinsic graphene and the highly doped graphene where 
the kink disappears. This finding is in agreement with 
the interpretation of S.Y. Zhou et al. [2ll o f the gap 
opening in graphene due to the substrate [22j and not 
to manybody effects. In fact even a small doping as the 
one present in the experiments of S.Y. Zhou et al. makes 
the kink too small to affect the ARPES results. Notice 
however that using a pertubartive solution for the Dyson 
equation, see Eq. [TJ we are not able to describe the plas- 
maron solution found by ref. [Til . [23j | , so we cannot ex- 
clude many-body effects in the ARPES results. Moreover 
the screeening is evaluated only in plasmon-pole approx- 
imation so effects due to presence of additional poles in 

are not considered (see ref. [2(j for a discussion). 
Then we studied the change in topology of the quasi- 
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particle band structure with the doping, in order to 
clarify recent experiments on doped graphene and on 
graphite that found a large trigonal distortion in the band 
structure 0]. In figdjwe report the vf along the direction 
K M and KT&s function of the doping in LDA and GW. 
As expected the GW results slowly go closer and closer 
to the LDA ones for increasing doping. Notice however 
that the change of the vf in the two directions as func- 
tion of the doping is different. In the inset wc report 
the ratio of the Fermi velocity along the line KM and 
K r. This value can be used as measure of the Dirac 
cone anisotropy. We found that the anisotropy is larger 
in GW and increases with the doping. A similar result 
was found also by Roldan et al. [24[ using tight-binding 
model, even if the increasing of the anisotropy they found 
is larger than ours. Finally we want to comment that in 
order to understand the size of trigonal wrapping found 
by ref. it is necessary to consider the hybridization of 
carbon tt bands with the dopant ones, that causes a dis- 
tortion of the band structure close to the Fermi energy, 
see for instance ref. 1251. 
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FIG. 3: Band gap of graphene at the M point in GW and 
LDA as function of the doping , divided for the LDA gap 
at zero doping. We report the calculations for two different 
distances of the graphene slab. 
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FIG. 4: Fermi velocity vf as function of the doping , in the 
inset we show the ratio of the Fermi velocity along the line 
If rand KM 



we found the presence of a kink in the graphene band 
structure, but it rapidly disappears with doping, at least 
in GqWq approximation using a plasmon-pole model for 
Wo- This finding is in agreement with the experimen- 
tal results and their interpretation by Zhou et al. [2lj | . 
Finally wc report the slope of the Fermi velocity in the 
directions KM and .KTas function of the doping, and 
found that the trigonal wrapping is enhanced by the GW 
corrections and it increases with the doping. 
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V. CONCLUSION 

We have studied the quasi-particlc band structure of 
graphene as function of the doping. The bands dispersion 
is obtained by the real part of the self-energy calculated 
in GqWq approximation, that is generally considered ac- 
curate for a broad range of materials, including graphene. 
We found a strong renormalization for intrinsic graphene, 
that rapidly decrease as function of the doping. This re- 
sult can be used as new reference for future experiments 
on doped graphene and intercalate graphite because it 
provides for the first time an estimation of the Fermi ve- 
locity change with the doping. Close to the Dirac point 



APPENDIX A: CONVERGENCE OF THE 
CALCULATIONS 

In a typical GW calculation different convergence pa- 
rameters have to be carefully considered, here we want 
stress the relevance of the convergence in the k-point grid 
and the distance between the graphene sheets. The major 
part of the GW calculations are performed using a su- 
percell approach, where in order to simulate an isolated 
system a large distance d between the slab is used. It is 
well known that the GW corrections for a 2D systems 
are larger because the screening of Coulomb potential is 
weaker making the electron-electron interaction stronger, 
but in a slab geometry this is true only in the limit of 
infinite distance between the slabs. Here wc report the 
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change in the GW corrections at M point as function of 
the distance fig. [5j However notice that in the experi- 
ments the screening by the substrate will never allow for 
the observation of a completely single sheet, so the direct 
comparison with experiments remains more difficult. As 



0.85 




0.45 u 1 ' 1 1 ' 1 ' u 

10 15 20 25 30 35 40 45 50 

distance between the slabs (a.u.) 

FIG. 5: Converge of GW correction to the gap at M as func- 
tion of the sheets distance. 

one can see from the fig. [5] in order to converge the calcu- 
lations with an accuracy of 0.05 eV a distance larger than 
50 a.u. is necessary, however in this paper we performed 
all calculations using a distance of 20 a.u. mainly for two 



reasons: first the behavior of the GW corrections with 
doping remains unaltered; second calculations with dis- 
tances equal or larger than 50 a.u. are computationally 
unfeasible using plane- waves. Anyway even if we cannot 
provide the full band structure with large distance we 
used a fix grid of k-points to extrapolate the value of the 
Fermi velocity for different sheet distances to the infinite 
one vp° . In table |T] we report the final result that was 
obtained fitting the vp with vp(d) = vp° + A/d, where 
d is the distance between the sheets and A and vp° are 
fitting parameters. 

Notice that as the distance between the sheets is in- 
creased the number of bands necessary to achieve the con- 
vergence in the screening W and in the expansion of the 
Green Function G rapidly increase as in a simple quan- 
tum well where the number of bands below a give energy 
threshold increases with the system size. Therefore we 
maintained a fix cutoff in energy for all the calculations 
and increase the number of bands when necessary. 
Another important parameter is the k-point grid. Here 
we report some convergences tests, see fig. [6j for the 
gap correction at 2Was function of the k-point sampling. 
We want underline two important aspects: first even if 
the convergence of the exchange part of the self-energy 
can be very slow fig. [6] panel 6), the combination of ex- 
change and correlation part converges quite rapidly with 
the k-points; second for doped graphene, the convergence 
becomes harder with increasing doping because a larger 
Fermi surface has to be sampled. 
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FIG. 6: a) Converge test of the GW correction to the gap at 
Appoint in graphene as function of the k-point grid plotted 
for different number of electrons. In x label we report the size 
of the grid, namely n stays for a grid n x n x 1. All results 
are converged within 0.01 eV. . b) convergence of the E x for 
the 7r band at the M point. 



